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■ Nucleon magnetic moments display a rich nonanalytic dependence on the quark mass in both 

quenched and full QCD. They provide a forum for a detailed examination of the connection between 
quenched and full QCD made possible through the formalism of finite-range regularised chiral ef- 
fective field theory. By defining meson-cloud and core contributions through the careful selection 
of a regularisation scale, one can correct the meson cloud of quenched QCD to make full QCD 
predictions. Whereas past success is based on unquenching the leading-order loop contributions, 
i r here we extend and test the formalism including next to leading-order (NLO) loop contributions. 

" We discuss the subtleties associated with working at NLO and illustrate the role of higher-order 

Cl^ corrections. 
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1 " I. INTRODUCTION 

>: 

(N. .„ 

, The study of the properties of hadrons continues to attract significant interest in the process of revealing and under- 
f^) ■ standing the essential mechanisms of QCD, the fundamental theory of the strong interactions. The non-perturbative 
lO \ properties of QCD and the difficulty of numerically simulating the theory at the light quark masses of Nature has 
provided a rich history of phenomenological models at both the quark and hadronic levels and an intense study of the 
quark-mass dependence of hadronic observables in effective field theory (EFT) . 

Based on the observation that all hadron properties show a slow, smooth variation with quark mass for pion masses 
above ^0.4 GeV, one can conclude that the nonanalytic contributions from pion loops are suppressed there 
An alternative regularization method, namely finite-range- regularization (FRR), resums the chiral expansion in a 
manner that suppresses loop contributions at large pion masses. Inspired by quark models 043 that account for the 
/\ \ finite-size of the nucleon as the source of the pion cloud, FRR EFT has been used to describe lattice data over a wide 
range of pion masses. 

FRR EFT was first applied in the extrapolation of the nucleon mass and magnetic moments 5- 7]. The remarkably 
improved convergence properties of the FRR expansion mean that lattice data at large pion masses can be described 
very well and the nucleon mass obtained at the physical pion mass compared favorably with the experimental value. 
Later, the FRR method was applied to extrapolate the vector meson mass, magnetic moments, magnetic form factors, 
strange form factors, charge radii, first moments of GPDs, etc. [§l-fl6j . The results are reasonable and reflect the manner 
in which FRR EFT characterizes the essential features of QCD at the hadronic level. 

The prevalence of the quenched approximation in the history of lattice QCD simulations provided an opportunity 
to explore the possible connection between quenched and full QCD data. Indeed quenched chiral perturbation theory 
was developed [l7l - [2l| to understand how the nonanalytic structure of quenched QCD differed from that of full QCD. 
It was the advent of FRR EFT that made it possible to define a pion-cloud contribution to hadronic observables and 
then proceed to correct the quenched cloud to that of full QCD |22| . 

The meson loop contributions are calculated in both quenched and full QCD in terms of the axial coupling constants. 
One then fits the quenched FRR EFT to lattice QCD to learn the low energy coefficients. This is done by fitting 
the coefficients of the residual series of terms analytic in the quark mass. With the assumption that the SU(3) axial 
coupling constants, F, D and C, do not differ significantly between quenched and full QCD, one can replace the 
quenched meson cloud contribution of FRR EFT with the full QCD cloud contribution. We note that in dimensional 
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regularization for example, the low-energy coefficients are composed of both residual series and loop contributions 
with no recourse to separating the origin of terms contributing to the total renormalised coefficient. On the contrary, 
in FRR the low energy coefficients of the residual analytic expansion provide the core contribution which is considered 
invariant in moving from quenched QCD to full QCD. 

Quenched and partially-quenched FRR chiral EFT has been used to stud y b aryon electromagnetic phenomena 
including charge radii, strange magnetic moments and strange form factors [12Hl5l ]. In the previous calculations, only 
the leading-order diagrams were included and unquenched. For example, power counting with Ma = Mjy for the 
magnetic form factor, only the leading nonanalytic terms proportional to log^n^) and m n were included. In this 
paper, we will include the next-to-leading-order (NLO) contributions. 

The paper is organized in the following way. In section II, we briefly introduce the relevant chiral Lagrangian. In 
section III, we study the nucleon magnetic moments using chiral perturbation theory with FRR at NLO. Numerical 
results and discussions are presented in section IV. Finally, section V provides a summary. 



II. CHIRAL LAGRANGIAN 



There are many papers which deal with heavy baryon chiral perturbation theory - for details see, for example, Refs. 
[23M25I ] . For completeness, we briefly introduce the formalism in this section. In heavy- baryon chiral perturbation 
theory, the lowest chiral Lagrangian for the baryon-meson interaction which will be used in the calculation of the 
nucleon magnetic moments, including the octet and decuplet baryons, is expressed as 

C v = iTrB v (v ■ V)B V + 2DTtB v S^{A^, B v } + 2FTrB v S v l [A ti , B v ] 

-if£(v ■ V)T Vfl + C{f^A^B v + B v AfjT^), (1) 

where 5„ is the covariant spin-operator defined as 



SH = l -j 5 a^v u . (2) 

Here, v" is the nucleon four velocity (in the rest frame, we have v" = (1,0)). D, F and C are the axial coupling 
constants. The chiral covariant derivative is written as D^B V = d^B v + \Vn,B v ], The pseudoscalar meson octet 
couples to the baryon field through the vector and axial vector combinations 

V, = \(Cd,C f + C + ^C), A, = i(C^C f - C + ^C), (3) 

where 

C = e** /f , / = 93 MeV. (4) 
The matrix of pseudoscalar fields <fi is expressed as 

7^ - + K+ \ 

(5) 

B v and T£ are the velocity dependent new fields which are related to the original baryon octet and decuplet fields B 
and T» by 

B v {x) = e lmN ^ v ^B{x), (6) 




T${x) = e mNl/v ^T' 1 {x). (7) 

In the chiral SU(3) limit, the octet baryons will have the same mass tub- In our calculation, we use the physical 
masses for the baryon octets and decuplets. The explicit form of the baryon octet is written as 
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For the baryon decuplet, the symmetric tensor carries three indices and is defined as 

T m =A++, T 112 = -^A+, T 122 = i=A , (9) 

222 — A , ^113 = -^=^ , J 123 — ' 
223 = ^^ , 133 _ ' 233 ~~ ^ ' i 333 — " ■ 

The octet, decuplet and octet-decuplet transition magnetic moment operators are needed in the one-loop calculation 
of nucleon magnetic form factors. The baryon octet magnetic Lagrangian is written as: 

C = ^- (fi D TrB v a^ {F+, B v } + ft F TvB v a^ [F+, B v ] ) , (10) 

where 

F+ = \ (C^F^QC + C^QC 1 ) • (11) 

Q is the charge matrix Q =diag{2/3, —1/3, — 1/3}. At the lowest order, the Lagrangian will generate the following 
nucleon magnetic moments: 

^ rcc = ^D+MF, »T = -\»d. (12) 

The decuplet magnetic moment operator is expressed as 

g _ _ 

c = -i—fJ-cqijkT^T^F^, (13) 

where qijk and qijkPc are the charge and magnetic moment of the decuplet baryon T^-fc. The transition magnetic 
operator is 

C = i^VtF» v U lk QlBi m SZT:- klm + ^ k Q[fl\ klm S» v Bl n \ . (14) 

In Ref. (26[, the authors used fi u , fid and fi s instead of the fic an d fix- For the particular choice, fi s = fid = — 5M«) 
one finds the following relationship: 

3 2 

(J>D = ^Mm, f*F = -f+D, fJ-C = fJ>D, = 4/i£). (15) 

In our numerical calculations, the above formulas are used and therefore all baryon magnetic moments are related to 
one parameter, fir>. 

In the heavy-baryon formalism, the propagators of the octet or decuplet baryon, j, are expressed as 

and I. XiN , ( 16 ) 



v-k- 5i N + ie v-k- 5i N + ie ' 

where is v^v" — g 11 " — (4/3)S'(, t 5^. S ab = — m a is the mass difference of between the two baryons. The 
propagator of meson j (j — ir, K, rj) is the usual free propagator, i.e. 

(17) 



k 2 - Mf + ie ' 



III. NUCLEON MAGNETIC MOMENTS 

In the heavy baryon formalism, the nucleon form factors are defined as: 

< B{p')\J^\B(p) >= u(p') \v,G E {Q 2 ) + ie ^P v ° S ^ Gm {Q 2 ) \ u(p), (18) 

I m N ' 
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Octet 

Decuplet 

7C,K,T| 

Y 

FIG. 1: Feynman diagrams for the nucleon magnetic moments. The last two diagrams, (g) and (h), only exist in the quenched 
case. 



where q = p' — p and Q 2 



According to the Lagrangian, the one-loop Feynman diagrams which contribute to 



the nucleon magnetic moments are plotted in Fig. 1. The intermediate baryons can be octets and decuplets. Diagrams 
(a) and (b) are for the leading order, while diagrams (c), (d), (e) and (f) enter at NLO. The last two diagrams exist 
only in the quenched case where the 7/ is degenerate with the pion and no if -meson loops contribute. 
The loop contribution to nucleon magnetic form factors at leading order is expressed as 



G 



■P(LO) 
M 



m N 
8^f 2 



NN tNN 
1 1 



(P) l7r 



oNA . 
' PlK( P yiK 



tNA , oWS t 



NT. 

K( P yiK 



oNA tNA 
' Pl7r(p) J l7r 



oNE* rNT' 
PlK(p) 1 lK 



(19) 



„n(LO) _ m N 



oNN T NN i oNT t NT i oNA T NA , oJVS* rATE* 
Pl7r(n)- , lir + Pi K(n) 1 IK + Pln-(n) 1 Itt + Pi K(n) 1 1 K 



The integration is expressed as 



ra0 



/ 



3u V 



d A k 



(20) 



(21) 



where 



Af = Ljj(t + t/2)u j (k 7 - Y/2)K(T + Y/2) + <5 Q/3 ) 



{u>j(t - Y/2) + 6 aP ){uj(it + t/2) + wjOt - Y/2)). 



(22) 



m 2 + fc 2 is the energy of the meson j. S a ^ is the mass difference between baryon a and (3. In our 



calculation we use finite-range regularization with u( k ) the ultra-violet regulator. This leading order contribution 
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TABLE I: Coefficients ft for the magnetic moments in full QCD case. 





oNN 

PlTV 


oNA 
PlK 


oATE 
PlK 


oiVA 

PlTT 


aiVE* 
PlK 


Proton 
Neutron 


(D + F) 2 
-(D + F) 2 


(D+3F) 2 
6 
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oJVEE* 
P5X 


oJVAE* 
PbK 


Proton 
Neutron 




4(D + F)C MT 


(D-F)(D+3F) 
6 

(D-F)(D+3F) 


5(D-F)C 
18 

(D-F)C 
18 


(D+ZF)C 
18 

(D+3F)C 
18 



has been studied in the previous paper which gives the leading analytic term to the magnetic moments. The first 
terms in Eqs. (|19[) and (1201) come from the tt meson cloud, while the last two terms correspond to the case where the 
intermediate baryons are decuplets. 

The NLO contribution to the form factors is expressed as 



G 



•p(NLO) 
M 



1 



48^/2 



oNN T NN , oiVS WV£ , dJVA tNA , aJVAS tAAE , oAA T NN , nJVA r AA , oiVE* r AE* 
P27r(p) i 2ir + P2K(p) 1 2K + l 3 2K(p) 1 2K + P5K(p) I 5K "+~ P2r;(p) J 2i) +P27r(p) J 27r + P2A(p)- , 2A' 



, nJVi rAA _,_ oASE* rJVES* , offAE* rJVAE* 
" i "P37i-(p) j 3tt T PhK(p) 1 bK ' PbK(p) 1 5K 



(23) 



G 



n(NLO) 
M 



where 



1 



jN'S* rAE* 



48tt3/2 

P37r(n) J 3ir 



oiViV TiViV . oiVZ^ riV2j , oiV A WV A . oiV A2j riV A2J . oiViV WV iV . oiVZl WV Zl . oiV 2j rj' 

P27r(n) i 2 7 r +P2A'(n) J 2A" +P2A(n)- , 2AT +P5A"(n) J 5A + P2r ) (n) J 2? ? T" P27r(n) J 2tt T P2Ar(n)- , 2if 

(24) 



oASE* rJVSS* 



,AA£* r JVAE* 
5A(n)- t 5A - "I" P5K(n) 1 5K 



1% = I d 3 k 



k 2 u 2 (k) 



(25) 



l 5j - 



k 2 u\k) 



Wj -(fc)( Wj (A) + <y^)( Wj (A) + <y«T))' 



(26) 



r"/3 
'3j 



= / 



k 2 u 2 (k) 



(27) 



All the coefficients /3 in front of the integrals are shown in Table I for full QCD. The coefficients of the leading-order 
contributions are functions of the coupling constants D, F and C. The coefficients of the NLO contribution are 
associated with the tree level baryon magnetic moments. 

The magnetic moment is defined as [i = Gm(Q 2 — 0). The total nucleon magnetic moments can be written as 



AV«) = ag + almi + a p A ml + (Z - l) M * ree + G 



p(LO) (Q 2 =0) + G^ NLO) (Q 2 = 0), 



M 



(28) 



fi n (ml ) = a£ + a>2 + ajm 4 . + (Z - 1 R ree 

where the wave function renormalization can be calculated as 

f 



G 



i(LO) 
M 



(Q 2 = 0) + G 



n(NLO) 
M 



(0 2 = o), 



i 



48^3/2 



NN T NN 
l 2j 



NA T NA 



NA T NA 
2j 



OiVA W' 

PK 1 2 



AS r AS 



PK 1 2j 



P 



NN 7-AS* 
K Hj 



p l ; 



NN T NN 
l 2j 



(29) 



(30) 



6 



TABLE II: Coefficients /3 and ft for the wave function renormalization in full QCD and quenched case. 



Full QCD 




ft"* 












9 j(D + Ff 


2C 2 


±(3F + D) a 




5 n2 
fi C 


i(3F-D) 2 
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oNN 


3JVJV 
Pdh 


oNN 
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-\D 2 - \F 2 + fDF 


Xff2 

9^ 


-§£> 2 - §F 2 + DF 


|Mo 2 (3F-_D) 2 


3(3F — D)(D — F) 





TABLE III: Coefficients ft for the magnetic moments in quenched case. 
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3JVA 
P2ir 



5NN 
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Proton 
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-^D 2 -±F 2 -^DF)y D 



3D A +3F*-2DF 



12 
-3D 2 



{^D + 3/J, F ) 



-3F*+2DF 



fl D 



5C 2 ,, 



(3F-D)(F-D) , . Q \ 
i ^ '(/XZ) + 3^f) 

2(3F-D)(D~F) 

3 MO 



The coefficients /3 in the wave function renormalization are listed in Table II. 

With the exception of Figs. [1] (a) and (b), the contributions of the diagrams in Fig. 1 are proportional to the 
tree- level moments, Mp^ ) expressed in Eq. (fl"2]) . In the quenched case [T(J, the logarithmic divergence of the magnetic 
moment encountered in the chiral limit makes it necessary to replace the leading order estimate M P r (^) with the 
renormalised moment, effectively incorporating physics associated with higher-order terms of the expansion. To 
provide a connection between the quenched and full QCD expansions, we make this replacement for the full QCD 
case as well. Therefore, the expression for nucleon magnetic moments can be written as 



p(n) 

pin) , pin) 2 , pin) 4 , pin) . / n -, \ , Mi!3 

Mp(n) = V + a 2 m 7r + V m 7r + Ma ' + {Z - l)Mp(n) + — ^p(n), 

Mp(ri) 



(31) 



where /ufi is the loop contribution from diagrams (a) and (b) in Fig. 1, while (Jr^,™ 1 is the contribution from (c), (d), 
(e) and (f) expressed in the previous formulas. The above formula can be rewritten as 



,p(™) ; c 



Mp(n) 



(™) , _P(") 



p(™)_„4 

04 y m„. 



f)}/(2 



p(n) 



M 



p(n) 



(32) 



Since the lattice data of the magnetic moment are obtained in the quenched approximation, we should fit the lattice 
data using quenched chiral perturbation theory. In the quenched case, only the pion loop makes a contribution. The 
coefficients in the quenched case are shown in table III. They can be obtained following the methodology of Ref. . 
Remember, in this case, we have two more diagrams, i.e. (g) and (h) in Fig. 1. 

The loop contribution to the nucleon magnetic moments at leading order in the quenched case is expressed as 



G 



P(LO) 
M 



m N 



oNN T NN , pNA T N, 



oNN tNN 
Pdh(p) 1 6ir 



(33) 



pn(LO) _ m N 
~ 87T3/2 



oNN T NN 1 qNA rJVA , oNN T NN 
Pl7i-(n) i l7r Pl7r(n) i l7r ~T Pdhfji) 1 ^ 



(34) 



where 



(35) 
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The NLO contribution can be written as 



G 



■p(NLO) 
M 



m N 

487T3/2 



qNN tNN 
j 2tt(p) 1 2tt 



NN tNN 



NA tNA 
27r(p) J 27T 



NN jNN 



(36) 



G 



n(NLO) 
M 



m N 
48tt3/2 



NN r NN 



riViV , o 

2it{n) 1 2-n: T P2rj(n) ± 2r) 



NN T NN 

1 9 



NA T NA 
2ir(n) 1 2TT 



■P. 



NN T NN 
h(n) 1 2n 



(37) 



In the quenched case, the wave function renormalization constant is obtained as Z 

1 



Z = 1 



48^3/2 



NN T NN 



1 2j Ptt *2j 



NA T NA qNNjNN 



NN T NN 



dh 6j 



Ph lf . 



2j 



(38) 



where the coefficients /3 are shown in table II. For the double hairpin diagram, Mq is the interaction strength. 
Similar to the full QCD case, the quenched magnetic moments of the nucleon are expressed as 



_ /„p(«) 
n) — y% 



pin) 



p(n)_4 



■/if 1 ( " ) }/(2 



~p(«) 

r.trcc 

M P („) 



(39) 



where /if^" is the loop contribution from diagrams (a) and (b) in Fig. 1 with quenched coefficients, while ^ s ^ ne 
contribution from the other diagrams. Because the simulation is on a lattice with length L in the spatial dimensions, 
the momentum integral is replaced by a discrete sum over the momentum, i.e., 



(40) 



where the momenta k x ,k y and k z are given by 2-KnjL and the infinite sum is regulated by the finite-range regulator. 
By fitting the quenched lattice data with Eq. (l39l) . one can get the parameters a^. The full QCD results are then 
obtained with Eq. ([52]) . 



IV. NUMERICAL RESULTS 



In the numerical calculations, the parameters are chosen as D = 0.76 and F = 0.50 (gA = D + F = 1.26). The 
coupling constant C is chosen to be —1.2 which is the same as Ref. [27|. The regulator u(k), may be chosen as a 
monopole, dipole or Gaussian function, since all have been shown to yield similar results [28j. In our calculations the 
dipole function is used: 

u{k) = (TTfcW' (41) 

with A = 0.8 GeV. 

The K- and 77-meson masses have relationships with the pion mass according to 

m K = \ m l + m K\phy- T^lphy, (42) 



m l = + m l Iphy - gm^lphy , (43) 

and enable a direct relationship between the meson dressings of the nucleon magnetic moments and the pion mass. 

We begin by considering nucleon magnetic moments from the CSSM Lattice Collaboration [3(3] • The leading order 
result of the proton magnetic moment versus is shown in Fig. 2. The solid line is for the finite- volume quenched- 
QCD fit and the dashed, dotted and dash-dotted lines are for the infinite- volume full QCD results of tree level, leading 
loop and sum of tree level and leading loop, respectively. One can see that quenched lattice results can be described 
very well in quenched chiral effective field theory. At the physical pion mass, the proton magnetic moment jl p is about 
2.25 which is significantly smaller than the experimental data. With the obtained fitting parameters Oj, the full QCD 
results are determined and illustrated in the figure. 
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FIG. 2: The proton magnetic moment as a function of the squared pion mass. The solid line illustrates the finite-volume 
quenched-QCD fit to the lattice results. The dashed, dotted and dash-dotted lines correspond to the infinite-volume full-QCD 
results at tree level, leading loop and sum of tree level and leading loop, respectively. 



In the quenched case, the loop contribution is small. While quenched-QCD coefficients of nonanalytic terms are 
typically smaller than in the full QCD case, the dominant effect here is that the momentum integration is replaced by 
the finite-volume sum. The loop contribution in full QCD gives the dominant curvature of the pion mass dependence 
of the proton magnetic moment. The proton magnetic moment in the full QCD case is significantly larger. At the 
physical pion mass, the proton magnetic moment in the full QCD case, /j, p , is approximately 2.95 /zjv which is similar 
to the experimental value, 2.79 /xjv- 

The leading order result for the neutron magnetic moment versus is shown in Fig. 3. Again, the finite- volume 
quenched-QCD lattice results are described very well by finite- volume finite-range regularised quenched chiral effective 
field theory. The curvature of the line is small. At the physical pion mass, the finite- volume quenched neutron magnetic 
moment is around —1.5. The associated full QCD results of tree level, leading loop and sum of tree level and leading 
loop are shown as well. Similar to the proton case, the loop contribution changes smoothly at large pion mass and 
drops quickly at small pion masses. The total value of neutron magnetic moment from our leading-order calculations 
is n n ~ —1.96 [an, similar to the physical value of —1.91 /in- 

The NLO result for the proton magnetic moment versus is shown in Fig. 4. The solid line is the finite- volume 
quenched-QCD result. The dashed, dotted, dash-dotted and dash-dot-dotted lines are for the infinite- volume full-QCD 
results of tree level, leading order loop, NLO loop and sum of tree and loop contribution, respectively. At NLO, the 
quenched lattice results continue to be described well by finite- volume quenched chiral effective field theory. However, 
at NLO, the approach to the chiral limit displays some downward curvature associated with the new wave-function 
renormalization contributions which appear only at NLO. The wave function renormalization constant Z decreases 
quickly at small pion mass. 

At the physical pion mass, the infinite volume tree-level contribution to the proton magnetic moment changes from 
2.30 to 1.10. The leading loop contribution at the physical pion mass is 0.65. The NLO loop contribution has a smaller 
curvature than the leading loop. It contributes 0.78 to the proton magnetic moment. The sum of tree level, leading 
loop and NLO loop contribution to the proton magnetic moment is 2.53 /in to be compared with the experimental 
value of 2.79 [In- 

The NLO result for neutron magnetic moment fi n versus is shown in Fig. 5. The meaning of the different types 
of lines are the same as for Fig. 4. Here the wave- function renormalization has a more subtle effect. As anticipated, the 
NLO loop contribution has a smaller curvature than the leading-order loop contribution. At the physical pion mass, 
the tree level, leading loop and NLO loop contribute to the neutron magnetic moment —0.52, —0.62 and —0.66 (In, 
respectively. The total neutron magnetic moment at NLO is — 1.80/xjv which remains close to the experimental value 
of —1.91 [in- 
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FIG. 3: The neutron magnetic moment as a function of the squared pion mass. The solid line illustrates the finite-volume 
quenched-QCD fit. The dashed, dotted and dash-dotted lines correspond to the infinite-volume full-QCD results of tree level, 
leading loop and sum of tree level and leading loop contribution, respectively. 



We should mention that when we calculate the NLO loop contribution the tadpole diagram is not included explicitly. 
That is, the tadpole contribution is handled by adjusting the parameters a! i to cii, so that: 

a' a + a' 2 ml + a!±m% + fi tad ~ a + a 2 m\ + a 4 m* , (44) 

where fitad is the tadpole contribution to the magnetic moments. In Ref. (ilj . the chiral extrapolation did explicitly 
include this tadpole diagram but the numerical results were almost the same if we refit the lattice data without 
this diagram. This means that, in practice, the new parameters, a i: can compensate the contribution of the tadpole 
diagram. In the present work we explored both the explicit and implicit inclusion of the tadpole term, with the 
numerical results clearly favoring the approximation where the fitting parameters on the right hand side of Eq. f|44[) 
are the same in the quenched and full QCD cases. 

Chiral symmetry can be realized in a number of ways, resulting in different forms for the effective Lagrangian. In 
Ref. [29l | , the authors applied two different Lagrangian densities incorporating chiral symmetry to the problem of pion- 
nucleon scattering, with the nucleon represented by an MIT bag. In one case the interaction was confined to the bag 
surface, where only a Yukawa type NNir interaction appeared. The other formulation involved a volume-interaction, 
where a contact term (four particle NNtttt term) is required. Their conclusion was that transforming from the surface 
interaction to the volume interaction amounts to summing the contribution from all excited intermediate states of the 
confined quarks. That is, the two formulations give equivalent results if excited intermediate states are included. One 
can also study the magnetic moments with the pseudoscalar nucleon-meson interaction where no tadpole diagram 
appears - c.f. Refs. 0-3- With this background we conclude that the tadpole contribution to the magnetic moments 
from the contact term corresponds to the contribution of diagram (c) in Fig. 1 summed over an infinite set of highly 
excited baryon states and phenomenologically this appears to be appropriately incorporated through Eq. (|44j) . 

V. SUMMARY 

We have extrapolated quenched lattice QCD results for nucleon magnetic moments extending into the chiral 
regime (3(| to the physical pion mass using finite- volume finite-range regularised chiral effective field theory. Here, the 
NLO contributions are included, with the numerical results showing that the quenched lattice results are described 
very well. By fitting quenched lattice data, the parameters can be obtained and using the dipole regulator param- 
eter of 0.8 GeV the full QCD results are predicted. The infinite-volume full QCD results obtained at the physical 
pion mass are in reasonable agreement with experiment at both leading order and NLO. Thus finite-range regularised 
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FIG. 4: The proton magnetic moment versus squared pion mass. The solid line illustrates the finite-volume quenched-QCD 
fit. The dashed, dotted, dash-dotted and dash-dot-dotted lines correspond to the infinite- volume full-QCD results of tree level, 
leading loop, NLO loop and sum of tree level and loop contribution, respectively. 



TABLE IV: The obtained coefficients a; and magnetic moments. 





ao 


0,2 


<J4 


-LO 


-NLO 


LO 


NLO 


Exp. 


Proton 
Neutron 


2.34 
-1.39 


-2.08 
1.17 


1.24 
-0.70 


2.35 
-1.45 


2.31 
-1.41 


2.95 
-1.96 


2.53 
-1.80 


2.79 
-1.91 



chiral effective field theory provides an effective formalism for extending chiral perturbation theory beyond the power- 
counting regime and connecting quenched QCD and full QCD in a quantitative manner. The parameters and results 
are summarised in Table IV. 

It is interesting how the NLO contributions come in a compensating fashion. While each NLO contribution displays 
significant curvature in the chiral regime, the net contribution is relatively smooth and otherwise easily compensated 
for by the residual series expansion. We expect that this qualitative behavior will continue as additional higher-order 
terms are introduced, as we are informed by the lattice QCD results displaying a smooth slowly-varying quark mass 
dependence. Indeed, it will be interesting to examine more physical quantities to gain a deeper understanding of the 
utility of hnite-range regularised chiral effective field theory. 
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